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TECHNICAL MEMORANDUM 1315 


DISPLACEMENT EFFECT OF THE LAMINAR BOUNDARY 
LAYER AND THE PRESSURE DRAG* 

By H. Gortler 

1. INTRODUCTION 


The following considerations refer to two-dimensional, not neces- 
sarily steady laminar movements of an infinitely extended fluid without 
free surfaces, for very small friction, about a body-*- immersed in this 
fluid; they deal with the mutual influence between the boundary layer 
developing in the proximity of the body and the outer flow which is 
practically free from friction. 

We denote by X and Y the coordinates of a Cartesian right-hand 
system of the plane fixed with respect to the body, by x and y the 
customary boundary layer coordinates defined in a zone near the wall 
(that is, an orthogonal right-hand system where x signifies the 
wall-arc length measured from a fixed contour point, y the wall distance 
measured positively from the wall toward the fluid) . 

Furthermore, R(x) is assumed to be the radius of curvature of 
the wall which we define in the manner that has become customary for 
boundary layer investigations, in contrast to the usual mathematical 
definition, as positive for walls which are convex with respect to 

O 

the fluid. Finally, u and v are assumed to be the velocity 
components of the boundary-layer flow, U and V -the corresponding 
components of the potential flow about the body in x or y direction, 

U and V the velocity components of the potential flow in X or Y 
direction. 


*"Verdrangungswirkung der laminaren Grenzschichten und 
Druckwider stand. " Ingenieur-Archiv, vol. 14, 19^4, pp. 286-305. 

^Tbe contour of the cylinder cross section in the flow plane 
should, if no other requirements are made explicitly, be of continuous 
curvature and free of multiple point singularities; it may go to infinity. 

2 ...... ... . . 

In order to have in the defined range of x, y a reversible 

one-to-one relation between X, Y and. x, y, the requirement must 
be made that there R + y £ 0 . This does not impose any further 
limitation, however, since according to presupposition in the boundary 
layer region (y/Rj <<1. 
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If the veil-known boundary layer omissions are permissible , the 
general hydrodynamic equations of motion are reduced to the following 
boundary-layer equations: 


Su 5u Su 

— + u — + v — 

c)t , <3x <3y 


2 

1 Sp + v d u 

p Bx oy^ 


(1.1a) 


c)u , <3v _ n 

§; 5? * 


(1.31) 


where t signifies the time, p the constant density, and v the 
constant kinematic viscosity of the fluid. Furthermore, p = p(x,t) 
is the pressure impressed on the boundary layer by the outer flow. 
Insofar as the outer flow noticeably deviates from the regular potential 
flow about the body, p(x,t) must be determined on the basis of special 
considerations or be taken from pressure measurements. In the present 
report, we consider only cases where the potential-theoretical pressure 
distribution in first approximation is sufficient for the boundary-layer 
calculation according to equations (l.l). In this boundary-layer theory 
approximation, the potential velocity prevailing at the "outer edge" of 
the boundary layer is replaced by the potential flow prevailing about 
the body itself, thus by U(x,0,t), and the pressure term in 
equation (l.la) is calculated from Euler's equation of motion with 
U(x,0,t) = UQ(x,t) as 


1 5p = + U SUg 

p 5x St 0 Sx 


( 1 . 2 ) 


The boundary condition required in addition to the adherence conditions 
u = v = 0 at the wall y = 0 is that U assume asymptotically for 
increasing y^~Re the value3 TJ()(x,t). (Re is the Reynolds number 

formed with a characteristic length and a characteristic velocity. ) 

The described ommissions, correct in the limit for indefinitely 
increasing Reynolds numbers, are only an approximation for the large 
but still finite Re-values for which laminar flow is plainly still 
possible. Particularly the connection between the boundary-layer flow 
thus calculated and the outer potential flow, assumed to be undisturbed, 

^The imaginary line, denoted above, in the customary manner of 
expression, as the "outer edge" of the boundary layer, runs at the 
distance from the wall y = S(x,t) ("boundary-layer thickness") at 
which u practically attains the value UQ(x,t) for the accuracy 

requirements in each case. 
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does not satisfy continuity. Thus the question arises: What actual 

course does the outer, practically frictionless flow which is modified 
by the boundary layer take? This problem and the related one regarding 
the calculation of the pressure drag will be discussed below. The 
considerations will be applied to an example in all details. 


2. DISPLACEMENT EFFECT OF THE BOUNDARY LAYER ON THE OUTER FLOW 


We assume a body which, from a state of rest, is set in motion 
relative to the surrounding fluid. In this case, the boundary layer 
gradually developing from the start of the movement will be very thin 
at first, so that for sufficiently small times of movement the regular 
potential flow about the body itself, as outer flow, may be taken as a 
basis for the boundary layer calculation in good approximation. The 
well-known investigation concerning the origin of the boundary layers 

(H. Blasius^), and later extensions^; ^ are based on this assumption. 
However, with increasing times t the boundary layer increases, for a 
finite Reynolds number, up to a noticeable thickness; in general, it 
will influence the outer flow by mass displacement away from the wall 
to a corresponding extent. For a corresponding improvement of the 
customary boundary-layer theory approximation sketched in section 1, 
one must therefore consider the mutual influence of the increasing 
boundary layer and the outer-flow. It suggests itself to attain such 
an improvement in first approximation by using at first the displacement 
effect on the basis of the original boundary-layer calculation, and 
hence determining a correspondingly corrected course of the outer flow. 
In a second step one would then have to calculate the boundary-layer 
flow anew, taking the above new outer flow as a basis. However, it will 
generally not be possible to carry out this improved boundary-layer 
calculation, using the new outer pressure distribution, according to the 
old scheme (equations (l.l) with the boundary conditions formulated 
there in the text), since the curvature effects, so far neglected in 
the boundary-layer equations, generally contribute amounts which can no 
longer be neglected within the scope of such an improved boundary-layer 
calculation. We shall return to this later on (section 4). 


^lasius, H.: Grenzschichten in Flussigkeiten mit kleiner Reibung 

(The Boundary Layers in Fluids with Little Friction) . Dissertation 
Gottingen, Zeitschrift fur Mathematik und Physik 56 , 1908, p. 1. 
(Available as NACA TM 1256. ) 

^Boltze, E.: Grenzschichten an Rotationskorpem in Flussigkeiten 

mit kleiner Reibung. (Boundary Layers on Bodies of Revolution in Fluids 
with Little Friction) . Diss. Gottingen, 1908. 

^Goldstein, S., and Rosenhead, L. : Proc. Cambridge Phil. Soc. 32, 

1936, p. 392. 
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The improved flow calculation aspired to is of decisive importance 
for the numerical determination of the pressure drag7. This becomes 
particularly obvious, if one visualizes a body of finite dimension 
transferred from a state of rest into a state of rectilinear -uniform 
movement relative to the surrounding fluid of infinite extension. For 
this terminal state the potential-theoretical pressure distribution, so 
far, in first approximation, taken as a basis for the boundary-layer 
calculation, does not yield any pressure drag (d'Alembert's paradox). 
Thus, small as the changes in pressure gradient at the body caused here 
by the boundary layer may be relative to the potential-theoretical 
pressure gradient, they alone constitute the pressure drag. 

We want to emphasize here once more that, as we discussed before 
in the introduction, we do not include in our considerations flows 
where, for instance, by pronounced separation of the boundary layer, 
the outer frictionless flow is considerably transformed compared to the 
potential flow about the body. Of course, precisely such processes are 
of foremost importance for the originating pressure drag in many flow 
problems. Such processes cannot be included within the scope of the 
following considerations which are based on the boundary -layer theo- 
retical approximation. However, we should like to refer here to a 
related report by M. Schwabe^ where the pressure distribution after 
completed boundary-layer separation is determined according to. an 
empirical formulation for the example of the circular cylinder set into 
rectilinear-uniform motion from a state of rest. The space taken up by 
the pair of vortices developing on the back of the cylinder after the 
separation is determined by observation and then simulated by calculation 
by superposition of a suitable time-dependent source-sink flow on the 
ordinary potential flow about the cylinder. 0he then obtains outside 
of this space a streamline pattern which corresponds well to actual 
conditions, and a pressure drag caused by the nonsteady acceleration 
fields. 

Furthermore, we want to point out that J. Pretsch^ developed an 
approximative method for theoretical determination of the pressure drag 


^"pressure drag" = the component in direction of the movement of 
the resultant pressure force on the entire body surface, also called 
(less correctly) "form drag”. Pressure drag + friction drag = total 
drag (presupposing non-existence of free surfaces). 

8schwabe, M.: Uber Druckermittlung in der nichtstationaren ebenen 

Stromung (On the Determination of Pressure in a Nonsteady Two-Dimensional 
Flow). Diss. Gottingen. Ing. -Archiv 6, 1935^ p. 3^. 

^Pretsch, J.: Zur theoretischen Berechnung des Prof ilwider standee 
(Concerning Theoretical Calculation of the Profile Drag). Diss. Gottingen, 
Jahrbuch 1938 der deutschen Luftfahrt fora chung, p. I 60. (Available as 
NACA TM 1009. ) 
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of profiles in a steady flow by talcing the displacement effect of the 
boundary layer into account. In a procedure similar to ours he pre- 
supposes that no separation of the boundary layer, or only a slight one, 
takes place; in the steady case this signifies, however, a limitation 
to slender profiles at a small angle of attack so that here the total 
profile drag is due in the greatest part to surface friction. Our 
consideration, related in the basic idea, refers at the outset to non- 
steady flows as well and aims chiefly at the problem of the drag origin 
in movements of cylindrical bodies of arbitrary profile from a state of 
rest for small times after the start of motion, naturally without being 
limited to these movements. Also, it follows a different method. 

Whereas Pretsch makes the additional assumption that the friction losses 
in the wake behind the body may be neglected, probably satisfied in 
good approximation in view of his presupposition quoted above ("Bodies 
of Small Drag"), and then is able to determine, within this scope, the 
pressure drag from momentum considerations in a simple and general 
manner without actually having to calculate the corrected outer potential 
flow in every single case, we abstain from this or a similar simplifying 
assumption because we set ourselves a problem of a different type. 


X 


\ In the following considerations we shall first deal with the first 

\ step of problem formulation, namely, the correction of the outer poten- 
) tial flow by consideration of the displacement effect (resulting from 
the boundary-layer flow calculated in first approximation), and shall 
discuss a few questions arising in case of nonsteady conditions. After 
,^that the second step which leads to the calculation of the pressure drag 
will be treated in general and carried out numerically on an example. 





3. CORRECTION OF TEE OUTER POTENTIAL FLOW 


First, one has to find an appropriate measure for the displacement 
effect of the boundary layer on the outer flow. For that purpose the 
well-known boundary-layer theoretical length presents itself which is 
called "displacement thickness" and denoted by 5*. 

The formulation of the boundary condition for u(x,y,t) for 
indefinitely increasing y ^Re was based on the conception that the 
potential theoretical velocity distribution for the large (though still 
finite) values of the Reynolds numbers of interest may be regarded 
practically as constant in first approximation for the very small thick- 
ness of the boundary layer and may, therefore, be put equal to U(x,0,t). 
The definition of the displacement thickness 6* is likewise based on 
the conception of this streamline approximation |u(x,y,t) = U^(x,t) 

and thus for reasons of continuity V(x,y,t) = -ydUg /<3x] . 
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The quantity S* is explained as follows: The fluid volume trans- 

ported at a fixed point x at a fixed time t between the wall y = 0 
and the outer edge y = &(x,t) of the boundary layer in unit t ime 

r d 

equals J u dy. The potential flow about the body would transport, if 

the streamline approximation described above were taken for a basis, the 
volume U Q (x,t)5 through the- same cross-section in unit time. The 

difference between these two quantities, that is, the loss in rate of 
flow per unit time caused by the viscosity effect, leads by virtue of 


U 0 & " f Q U ^ = u o 5 * 


to the definition of the length 5* 


6*(x,t) 



(3.1) 


[Compare figure 1; the hatched areas are equal. Of course, any other 
length y = yj_ > 5 may be selected instead of the upper limit 5 of 

the integral in equation (3.1), as long as the velocity U(x,y,t) 
in 0 < y ^ yi is replaced by U 0 (x,t)]J According to definition, 

&* is, therefore, a measure of the displacement of the streamlines of 
the outer potential flow away from the body on the basis of the reduction 
(caused by the friction layer near the body) in the quantity of fluid 
flowing by the point x at the time t. (Therein a streamline in its 
identity for all times is prescribed by the fact that it, together with 
the wall y = 0, includes a stream tube of temporally constant through 
flow. ) 

In order to obtain, instead of the potential flow about the pre- 
scribed body, a corrected outer potential flow which takes the dis- 
placement effect of the boundary layer into account, we visualize the 
following model flow: Outside of a line y = S*(x,t) a potential flow 

with y = 5* as streamline is assumed to flow. Within 0 = y ^ 
one assumes water at rest relative to the body. In this model flow the 
same quantity of fluid is to flow past the body per unit time as in the 
viscous flov. Equation (3.1) then indicates the value of 5*(x,t) in 
first approximation for sufficiently high Re values. 
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As to the line y = 5* which is to he a streamline of our substitute 
flow, it must be taken into consideration that for the general nonsteady 
case the course of this line varies with time. The line y = 5* then 
is a movable dividing line between potential flow and water at rest. 

In contrast to the dynamic dividing lines, it therefore consists in 
general not permanently of the same fluid particles, but represents, 
for reasons of continuity, a permeable line. In mathematical formulation 
the boundary condition to be stipulated expresses that the normal 
component of the potential-flow velocity along the dividing line 
y = S*(x,t) should vanish every moment. Then y = 5* is a streamline, 

A simple example which we shall investigate more thoroughly later 
(section 5) will serve to clear up these conditions. We assume that a 
circular cylinder of the radius R is set at the time t = 0, from a 
state of rest, relative to the surrounding infinitely extended fluid, 
into a rectilinear and uniformly accelerated motion perpendicular to 
its axis. The frictionless flow relative to the body is then given by 
the velocity potential 


4 > 


0 



■0 


where r and ■d signify polar coordinates in the flow plane, referred 
to the center of the circle, r=R + y, d = it - x/R, and b denotes 
a constant acceleration. This potential flow yields the pressure 
variation on r = R required for the usual boundary-layer calculation. 
For very small times ( small compared to the time t = t^ of the start 

of separation) this calculation results in a displacement thickness 
5* increasing proportionally to 'Jvt. For these very small times after 
the start of the motion the potential of our improved outer flow 
therefore reads 


<X> = bt 


r + 


(B 


+ c 


yvfc ) 2 


cos t3 


(t > 0, r > R + c yvt; 
5* = c Vvt, c = const.) 
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The line y = c *Jvt is a streamline. While it travels out into the 
fluid starting from t = 0, the streamline pattern of the entire outer 
flow correspondingly varies continuously! 0 . 

In order to clarify the conception we shall add a few remarks 
regarding the mass displacement of the boundary layers. The streamline 
displacement thickness &* does by no means always represent at the 
same time a measure for the mass displacement. This becomes Immediately 
clear from the following simple example. An unlimited plane wall which 
up to the time t = 0 is at rest relative to the surrounding infinitely 
extended fluid is assumed to be moved in itself, according to an arbitrary 
acceleration law, starting from t = 0. A certain velocity profile 
(which is known to be easily defined) develops, and one obtains a dis- 
placement thickness S* different from zero. However, a mass displace- 
ment away from the wall does not take place due to v = 0 (u = u(y,t)] . 

We could have made clear the difference in the two displacement phenomena 
in a perfectly analogous manner on the case of the circular cylinder 
set rotating about its axis from a state of rest; thus we should have 
spared ourselves the fiction of an infinitely extended body. Here, just 
as in the above limiting case of the plane wall, a boundary layer 
develops; a mass displacement away from the cylinder in radial direction, 
however, cannot take place due to reasons of continuity and symmetry. 

If one relinquishes this symmetry by selecting for instance instead of 
the circular cylinder a cylinder slightly wavy in comparison, the mass 
displacement normal to the contour will be, in general, different from 
zero; however, the length 5* does not present a measure for it. 


10 A 
We shall give here for comparison the velocity potential <5 of 

the corresponding flow about an expanding impermeable circular cylinder. 

If its contour is prescribed by r = R(t), one has 


3 > 


bt 


r + 


R 2 (t) 


cos t3 


+ R 


f vdR(t) 



The additively appearing source yields that additional fluid in the 
outer space r = R(t) which the impermeable circular cylinder itself 
displaces while expanding 


/ _ dR (t) 
^ dr '' l r=R(t) dt 
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Looking for a measure for the mass displacement, one may start from 
the relation following from equation (3.1) by differentiation with respect 
to x and consideration of the continuity equation 

(U 0 B*) = v(x,8,t) + S (x,t) (3.2) 


(Again any y = y ^ >6 may be substituted for 5) . The right side 

indicates the excess of the v-velocity at the edge of the boundary 
layer compared to the V-velocity - 8c)Uq/<3x of the corresponding 

frictionless flow prevailing there. Thus the entire volume displaced 
by friction effect between two points x Q and x^ per unit time is 

at any rate always prescribed by 



- V) dx = 



We now visualize, as in the definition of &*, a model flow. We assume 
a frictionless potential flow, with the quantity 8^ generally different 

from 8*, to be prevailing for y = S^(x,t); the length 8^ is assumed 
to be fixed so that precisely the fluid which actually is displaced by 
the friction effect would flow through 0 = y = 5^ with the 
velocity Uq. Inside 0 = y = 8^ we assume for our model flow water 
at rest relative to the body. Then there follows from equation (3.2) 

£ ( u o s .) = s (D ° 8 * ) • (3 - 3) 


Hence there results by integration for the "mass displacement 
thickness" 8^ the statement 

V 5 * - & *) = f (t), (3.4) 

that is, this difference is only dependent on t. 

Hence one can recognize: If one were to proceed in the determination 

of the first correction of the outer potential flow described above. 
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and therewith of the pressure field correction for large Re-values, 
in such a manner that one would make not y = 8* hut y = 8 # for all 

times the streamline of the substitute flow, the result would remain 
unchanged, according to equation (3.4). Both models differ at any 
moment only in the numbering of the streamlines. 

As to the relation between 6* and 8^, one may state: If one 

assumes that at the respective time t the stream tube which includes 

at the point x the region 0 = y = y ^ with y^_ = & followed up 
upstream finally blends completely in a region where the flow is 
frictionless (case of approach flow), 8* = & . If, however, this 

presupposition is not fulfilled, as for instance in the examples with 
vanishing 8^ selected above, the conclusion drawn above is not valid, 
either. In every flow of this type the streamtube just considered is 
either bounded by the wall in its entire course upstream, or it leaves 
the wall in stretches as when a separation of the boundary layer material 
from the body with subsequent readherence takes placed- upstream 
(case of longitudinal flow.) The value of 8^ at the time t can be 

given for all x if one is able to give, in addition to the variation 
of S*(x,t), known according to equation (3.1) , the value of 8 # at a 

single point x = XqJ 12 for equation (3.4) then yields only the 
difference statement 


V* 


x 

*0 


u 0 6 * 


JO 


Tit is always presupposed that the flow in the x-interval of 
interest of the body contour and outside of the boundary layer next to 
the body is frictionless. Thus for instance the case where the body 
gets into its own wake during the motion is left out of consideration. 

1 P 

It must be noted, though, that cases exist where it is impossible 
on principle to give such a value. Let us visualize for Instance the 
boundary-layer flow originating when an unlimited wall, deviating from 
a mean plane by surface waviness, is moved in its mean plane out of a 
state of rest. If no special conditions prevail, the expression for the 
total volume displaced at the time t up to a point x in unit time 


/ 

— CO 


/ X du n . 7 

v(x,y x ,t) + y 1 (x, t ) 


dx 


will not even be unique. 
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Therewith we shall close this little digression concerning mass 
displacement. Following we shall make use only of the streamline 
displacement thickness 6*. 


4. CALCULATION OF THE PRESSURE BRAG IN FIRST APPROXIMATION 


With the improvement in the calculation of the outer potential flow 
a more accurate specification of the pressure distribution in the outer 
frictionless flow becomes possible. However, we must warn here against 
the following fallacy: In general, one cannot use the resulting pressure 

distribution at the edge of the boundary layer as impressed 
'-pressure p(x,t) for an improved boundary-layer calculation according 
to the equations (l.l) and (1.2), for the obtained correction of the 
outer pressure gradient is, in general, of the same order of magnitude 
as the terms which have been neglected in the equations of motion. To 
obtain an improved boundary- layer calculation it would thus be necessary 
to take corresponding further terms into consideration in the equations 
of-motion as well. 


The general Navier-Stokes equations of motion and the continuity 
equation read in our curvilinear coordinates x,y in full strictness^ 


du , R du chi uv 
— + u — + v — + 


R dp 


dt R + y dx dy R + y 


p R + y dx 


du 


+ 


+ v 


R 


-\2 
d u 

dy 2 


d 2 u 


R + y dy (R+y ) 2 dx 2 


u 


2R 


(R + y ) 2 (R + y)^ dx 


d v k 


R 


dR 


(R + y)^ dx 


v + 


Ry 


dR du 


(R + y)3 dx dx 




(4.1a) 


13 

Compare, for instance: Tollmien, W.: Grenzschichttheorie 

(Boundary-Layer. Theory) . Handbuch der Experimentalphysik (Manual of 
Experimental Physics) Vol. 4, Part I, p. 248, Leipzig 1931. 
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St St 
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+ 'v 


2R Su 


S^r 

Sy 2 (R + y) 2 Sx 


_ Sfr 

R + y Sy (R + y) 2 Sx 2 


_1 St + R^ 


R dR 


(R + y) 2 (R + y) 3 dx 


u + 


(4.2a) 


Ry dR St 


(R + y) 3 dx Sx_ 


E + + , 0. 


R + y Sx Sy R + y 


(4.3a) 


Therein R(x) denotes the radius of currature of the wall contour y = 0 
(contrary to the mathematical definition positive on walls conTex with 
respect to the flow). 


If one considers in these equations, on the basis of the customary 
estimates founded on the physical picture, only the terms of highest 
order for Tery large Reynolds numbers Re = UL/v (U is a characteristic 
Telocity, L a characteristic length), one arrives, in the known 
manner, at the boundary- layer equations (l.l)l4. if one considers in 
improTed approximation also the terms of the order 0(S/L) compared 
with 0(l), there results 


2 

cta+u — + v— +— — S u _ y ^ Su ut | y tip y Su 
St Sx Sy p Sx Sy 2 R Sx R pR Sx R Sy 


1 Sp u 2 

P Sy R ’ 


(4.2) 


14, 


'Compare, for instance. 


W. Tollmien, footnote 13. 
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(^. 3 ) 

bx by R Sx R 

On the left side are the terms from equation (l.l), on the right the 
newly added terms. 

The correction of the outer pressure we obtained above is, in 
general, of this same order 0(&/L). Thus, if one wants to perform with 
this improved outer flow an improved boundary-layer calculation (second 
approximation for large Re), one can in general no longer neglect the 
right sides and can, therefore, no longer calculate with a 
pressure p = p(x,t) impressed on the boundary layer. (Compare addi- 
ctional remark at the end of this section. ) 

For calculation, of the pressure drag in first approximation for 
large Re,' the problem that interests us here, the solution is simpler. 

We now know, according to the expositions of section 3, the outer pressure 
field, far remote from the body, sufficiently accurately; it only remains 
for us to continue this pressure field up to the body surface by deter- 
mining the pressure gradient in y-direction through the boundary layer. 
Equation (4.2) serves this purpose: It yields dp/Sy sufficiently 

accurately for this approximation, if we substitute in it on the right 
side u from the first boundary-layer approximation. 

If y = B(x, t) denotes the "outer edge" of the boundary layer, 
and pjx,B(x,t) ,tj is the pressure distribution of the improved outer 
potential flow along this line, one has for 0 = y = 5 in this 
approximat ion 


P(x,y,t) = p(x,B,t) - 


u dy. 


(4.4) 


The pressure drag of the unit length of the cylinder in the flow becomes 
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Wherein 9 signifies the angle between surface normal and main flow 
direction, and the integral is to he formed over the entire contour K 
of the cylinder cross section . ^ 


■'■^Additional remark at the time of proof correction: Regarding 

the problem of an improved calculation of the boundary layer flow 
(second approximation for moderately large Reynolds numbers), not 
followed up further in the present report, the following calculation 
procedure seems to me to be promising: 

1. Calculation of the boundary layer in the customary first 
approximation for very large Re. 

2. Hence correction of the outer frictionless flow according to 
Ziff's method. 

3. Improved calculation of the pressure field p(x,y,t) in the 
boundary-layer zone according to equation (4.4). 

4. Again calculation of the velocity components u and v of the 
boundary layer in second approximation from equations (4.1) and (4.3), 
using the pressure field calculated above and replacing the right sides 
of equations (4.1) and (4.3) by the known expressions of first approxi- 
mation so that the newly added terms of the order 0(d/R) compared with 
1 appear in the calculation as prescribed functions. 
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5. EXAMPLE 16 


Following we shall consider as an example a nonsteady flow which 
originates if a "body is, from a state of rest, set into a rectilinear 
motion relative to the surrounding infinitely extended fluid, onward 
from t = 0. The relative velocity of the undisturbed approach flow 
with respect to the body visualized as being at rest is assumed to 


XQ The simplest, almost trivial example on which the method developed 
may be tested is the case of the plane steady stagnation point flow. 

Here the strict solution of the Navier-Stokes equations is known, and 
the flow near the wall calculated on the basis of the boundary-layer 
theory is known to agree with the exact solution. If one replaces the 
potential-theoretical stagnation point flow at the wall by the stag- 
nation point flow at the wall shifted by o*(= const.), one obtains 
as the corrected outer flow for y > 5 also full agreement with the 
exact solution. Since the outer pressure gradient parallel to the wall 
has remained unchanged in this correction and the wall is plane, an 
improvement of the boundary-layer calculation proves to be impossible 
as it has to be. Thus the exact solution has already been attained with 
this one step. 

A further example with plane body boundaries (for which the 
boundary-layer equations (l.l) therefore are valid except for terms of 
the order 0(S 2 /R 2 ) compared with l) is the case of longitudinal flow 
over a plate. According to Blasius, elsewhere, 6* = 1.73 ^/vx/u^ 

(Uoo free stream velocity, x distance from the leading edge of the 
plate) . With the improved outer pressure distribution calculated as 
suggested above, there results as the superimposed pressure for an 
improved boundary calculation 

p(x) = jj- + 1.35 Re(x)J 

(Re(x) = Uoox/v). One does obtain here a considerable pressure correc- 
tion since the new outer flow has a stagnation point, but one recognizes 
that the pressure correction has already dropped at Re(x) = 75 to 
1 percent of the stagnation pressure. However, in the proximity of the 
leading edge of the plate (small Re(x) - values) Blasius' boundary 
layer equation cannot be used. An improved calculation of the flow 
which in this range does not use the boundary layer approximations is 
still lacking. Thus for the time being the method of improvement for 
larger Re(x) suggested in this report cannot be utilized for this 
example. 
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"be U = C^t 11 (c Q = const., n = 0, t = d). Then the potential 

velocity (x,t) along the hody boimdaries can he represented, in the 

form 


U Q (x,t) 


0 for 
q(x)t n 



for t = 0 [ 


(5.1) 


We assume ty(x,y,t) to be the stream function of the boundary-layer 
flow defined by u = cW/dy, v = - d\|r/dx which is obtained when the 
pressure gradient to be calculated from equation (5.1) according to 
equation (1.2) is taken as a basis, on the strength of the boundary- 
layer equations (l.l), thus in the customary first approximation for 
very large Re-values. Generalizing the series developments set up 
by Blasius^ for the special cases n = 0 ( sudden transition from state 

of rest to motion at constant velocity) and n = 1 (uniform acceleration) 
one obtains the result 


lKx,y,t) = 2fvt qt n £ n ^ 0 (T)) + iq , t 2n+1 ^ n ^ 1 (Ti) + ...1, (5.2) 


n~f*l 

thus a series development in powers of t u (that is, in powers of 

the distance covered by the body) . The appearing coefficient functions 

£ (ti), £ (t|), ••• with T| = y/2 ^vt are universal functions of tj. 

n,0 n,l 

More details may be found in the appendix to this report. 

From equation (5-2) follows 


, n y | 1 1 2n _ i - lt r 1 , 

U = q ^n,0 + qq S n,l + 

V = - 2 V^jq't 11 ^ Q + (q ;2 + qq")t 2n+1 £ n ^ 1 + 



(5.3) 


Compare footnote 4 on page 3. 
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The boundary conditions are satisfied by virtue of 


~ ' "• =°V 


E ’ „(o) - 5' (0) = ... = 0, 

n,0 n,l 


} 


= 1 ’ ? n,l W ‘ •" = °» 


where there is always, here and later on 


t' = a t 

n, k d.T] n , k 


For the streamline displacement thickness 5* there results 


S*(x,t) = 2 /\/vt a n - q.' (x)t n+1 ^ n ^ 1 (oo) 


3 


with 


a 


n 


= lim p 

T) -» oo 



(A few numerical values are given in the appendix. ) 

We consider as an example the flow about a circular cylinder, 
the mode of notation of section 3 , the velocity potential of the 
frictionless flow about the circular cylinder is given by 


2 

0 n = U(t)(r + — ) cos 
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Hence the velocity at the periphery of the cylinder is calculated as 

U Q (x,t) = q(x)t n = 2U(t) sin -d f (5.8) 

and the resulting pressure distribution at the body surface is 

+ 2U 2 (t) sin 2 J , (5.9) 

where p^(t) is the pressure at the forward stagnation point. The 

potential- theoretical pressure drag of the length L of the cylinder 
at rest is therefore 


Po(x,t) = p d 


- P 2 


2R(cos i3 + l) — 
dt 




cos t 3 Rd$ = 2pjtR 2 L — 

dt 


(5.10) 


Thus one obtains, as is known, an increase in inertia by double the 
amount of the inert mass of the fluid displaced by the cylinder itself. 
For the improved calculation of the outer potential flow there results 


S*(x,t) = 2 *fvt 


a n 


+ 2 


i , (“) t: 

R n,l 


n+1 


cos -5 + 


(5.11) 


Following we limit ourselves to small times t for which the two first 
terms of the series development represent a sufficient approximation. 

In this approximation the line y = S*(x,t) practically represents a 
circle; that is to say, a circle with increasing radius r, the center 
of which, for t>0, does not coincide with that of the cylinder cross- 
section, but travels slowly downstream. The equation of such a circle 
reads formally 


r = b cos 3 + a 


(- 


4 sin 2 «) 2 
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(a = radius, b = displacement of the center), and if b « a 


r = a + b cos $ 



We identify 


a(t) = R + 2a n *Jvt, 


b(t) = 


4»s/vt 

R 


£ (oo)C 0 t n+1 
n,l 



( 5 . 12 ) 


If, temporarily, ’r,a’ t denote polar coordinates about the center 
of this circle in motion relative to the cylinder, the potential flow 
appertaining to the latter, the streamline of which is the line 
r = R + &*, that is, r = a(t), is given by the potential 


®i = ( u - HX 1 + §I) r cos * 


thus the potential flow here required in the polar coordinate system r,-9 
fixed relative to the body by 


” 1 + 7 


At first we introduce dimensionless quantities as follows: We 

choose R as the characteristic length, and the time T the body requires 
to cover the distance R, starting from the beginning of motion t = 0, 
as the characteristic time, thus 

1 

n+1 


(n 


1} F 


(r cos d - b) 


- 2br cos -a + b 


(5.13) 


T 


(5.14) 
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Hence results the characteristic 


velocity 


R 

T 



and the Reynolds number Re formed with this velocity and the length 
R become sl8 


Re 



(5.15) 


By making dimensionless the variable lengths x, y } r } a., b by 
dividing by R, the time t by dividing by T, the velocity U(t) 
by dividing by R/T, the potential by dividing by R 2 /T, and the 
pressure by dividing by pR2/T2 (we denote the dimensionless quantities 
by adding a wavy line ) , we obtain 


_1 1 

a = 1 + 20^ Re 2 1 
^ ~n+2 f 

b = 4(n + l)£ («>) Re 2 t 2 j 

n,l 


(5.12a) 


* = - (n + l)t n 1 + 


r 2 


+ 2 b r cos x + b‘ 


(x cos x + b ) . 


(5.13a) 


18- 


'If one forms the Reynolds number Re, with the length Wt 



one obtains 
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Furthermore, 


= 


2 i^vt 


I Rek-i 


and 


1 

= ^“x|2a n - 4(n + l)^ n 1 (~)t n+1 cos x~j. (5. Ha) 

Re^ L ' J 


We consider first the case n = 1 of uniform acceleration from a 
state of rest. Following, P.^ without more precise data will represent 

those additive portions of the pressure p^ = p^/(pR 2 /T^) which do not 

make any contributions to the pressure and are therefore not of interest 
to us. Furthermore, terms of the order 0(S 2 /R 2 ) compared with 1 are 
neglected. One obtains for the pressure calculation according to the 
Bernoulli equation at the distance from the wall y = 6/R 


5 

St 


^ = 4 (1 + 3a Re 2 t 5 ’) cos x + P + o( §E.''| 

£ \ 1 / 1 \ r 2 J 


1K 2 q/'b 2 ' 

2\3r / " 


- *4 . 16 £ t2 Bln^ + P„ + of ¥ 


2\? ax 


R 




Thus one obtains 


P]_^x, tj = 4(1 + 3<XjRe""^£^ cos x + l6 t^sin^x + P^ + 


(5.16) 
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According to equation (4.4) 


o, V) , S^X, |, t) - /“ 2 ay- 

0 


(4. 4a) 


In the small times of interest to us there is 


u = 4 £ sin x£' (tj) — l6 ^sin x cos x£ 1 . (q). 


1,0 


1,1 


therefore 


8 /r , 2 r s/B , 2 

-i ^ ^ 0^^ + 


u dy = - 16 t-^sin x 
0 ° 0 


128 t^sin 2 x cos x / t' t' dy - 

J 0 1,0 1,1 


■b/R 


P S/R 

256 t^sin^x cos 2 x / £’ 2 dy. 

J 0 1,1 ^ 

"(5.17) 

Because of lim £' (q) = 1, lim t,' (q) = 0 the two last integrals 

T)-> co 1,0 tj 00 1,1 

on the right are, if S/R is chosen so large that these asymptotic 
values are attained with the desirable accuracy, numbers independent of 
the variation of the boundary- layer edge S(x,t) with x. The first 
integral at the right, together with the corresponding second term at 
the right in equation ( 5 . 16 ), may also be combined into an expression 
independent of S(x,£), namely the expression 


16 t^sin^x / (l - £-! Q^)dy . 
J 0 3 
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Thus the variation of 5(x,^) over x does not play any role in 
judging what terms make contributions to the pressure drag which is as 
it should be. If one finally inserts 


_1 1 

dy - 2 Re 2 t 2 dq, 


all together one obtains, therefore, 

1 IS 


p(x, 0, £) = 4^1 


„ „ 2~2 , ^ 
+ 3ct^Re t ^fcos x 


•X 


-i 2 


+ 256 sin% cos x Re 2 t 2 / Q ^ + OCRe" 1 ). 

J 0 ’ 3 

(5.13) 


Thus the following- pressure drag of a circular cylinder of the length L 
for small times after start of the motion and in first order for large 
Re -values results 


W. 


D 


2priR 2 LCQ 


= 1 + 


1 JL 

^Rr 2 


(3^ + 


16 ^ 




(5.19a) 


2p tR LCq = is the potential-theoretical pressure drag. (Compare 

equation (5.10).) According to Blasius 'calculations 3ccq = 2 /tfr( =1.128), 
furthermore, according to the author's numerical evaluation 


Si.O 5 !,!* 1 - °- 09e ° k 


0 
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Thus the result in dimensional form is 



D r 


1 + 


Vvt / 2 , 

R \i[7 


0.392 



As friction drag 


(5.19) 




t) 


sin — cbc 
R 


there results from Blasius' 


calculation results 


19 


W R = 



2 

r yV* 


0.029 


R‘ 



(5.20) 


For sufficiently 

is, the friction drag 
according to the same 
the pressure drag. 


small times t(t « l) thus 


W R = W D “ \> that 


increases immediately after start of the motion 
law as the contribution of the skin friction to 


The total drag becomes 


W = U D + W R " V 



(5.21) 


In figure 2 we represented these results using the dimensionless 
"t = t/T (with T = ^2R/Cq for the present case n = 1), thus the 
relations 

tr 2 - = 1 + “ ^ #(1.128 + 1.569 t 4 ), (5.19-,) 

•'D n R X 


-*-9xn order to obtain the term with t , one must here include in 
equation (5.2) also the third term of the series. The necessary data 
may be found in Blasius’ report, footnote 4. 
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= 4^. t 5 (l.i28 - 0.1l6t 4 ) , (5.20 ) 

X R 1 


_H_ = 1 + ^( 2,257 + 1.^53^). (5.21-.) 

% R 


Our formulas can "be applied with good approximation only for very small 
times after start of the motion (solidly drawn parts of the curves), 
because we had, for the sake of simplicity, taken into consideration 
only a few terms of the developement (5.2); hut of course, with a little 
calculation expenditure they can he easily improved, at least so far 
that they are valid up to times shortly after setting-in of the separation. 

The value t = t^( = . 586 ) = 0.7 66 plotted in figure 2 indicates 

the time of the start of separation, in first approximation, at the 
rearward stagnation point according to Blasius. (Compare also appendix.) 

For an arbitrary integral n > 0 one obtains correspondingly as 
pressure drag with n(t) = Cgt 11 and Wj-^ = 2pitR 2 L du/dt 


W 


D 


Wn 


'0 




2n+2 


1 ^ S' 1 dr) 

0 n,0 n,l 


( 5 . 22 ) 


For the calculation of the coefficient of t 2n+2 sufficient numerical 

data concerning £ 1 (t)) are lacking so far. For this reason we limit 

n, 1 

ourselves in our numerical statements to the first term of the develop- 

1 

ment with time (term with t 2 ) and put the question whether the law 
W R = W D - W Eq obtained above for n = 1, for times immediately after 

start of the motion, is valid also for arbitrary n > 0. The friction 
drag has for these small times the value 


Wt, = W- 


D r 




R 


2n 


(5.23) 
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In the appendix the exact expressions for 
However, the relation 


i 



are derived. 


2(2n + l)an = ^,0 (0) 


2 2n (n j ) 2 2 

(2n) • tfn 


( 5 . 24 ) 


is valid; we produce a very simple proof for it in the appendix (one 
may confirm it also with the aid of the tables of the appendix we 
calculated for n = 1, 2, 3, and 4). Thus, for very small times after 
start of the motion, there applies indeed 


W-r = Wj-, - W 


D 


v D r 


viz : 


= W- 


Dr 


Vvt 2 2n ni(n - l) 
R ( 2 n)I \fn 


> 


(n > 0 ) . 


(5.25) 


(One has 2 2n n.'(n - 1)1/ (2a) l Vn = 1.1284 for n = 1; 0.7522 for n = 2; 
0.6018 for n = 3; O.5158 for n = 4, etc.) 


One may now consider more general laws of motion of the form 


r 


U(t) = 


0 


f(t) 


for t = 0 

for t = 0 (f( 0 ) = 0 ) 


and carry out corresponding calculations. If one presupposes that the 
function f(t) defined in t ^ 0 can he developed into a Taylor series 
aro un d t = 0 which converges for the small times 0 < t « t^ after 

start of the motion which are of interest to us, one may attain the 
result quite analogously with a series expression correspondingly 
generalized compared to equation (5.2). One can interpret the 

law U = Cgt n which has been valid so far as the first term of the 
development with time of such a general law of motion. Hence it follows 
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that for all these laws of motion the portion of the pressure drag 
caused by friction W-p - Wj^ increases immediately after the start of 

motion according to the same law ( 5 . 25 ) as the friction drag. This is 
a noteworthy quality of the circular cylinder. 


One may also include the case of a sudden start of motion in these 
considerations; it is true that one must then accept, corresponding to 
this degeneration of the form of motion at the time t = 0> infinite 
pressures and drags at the time t = 0. With U = Cq = const, for t = 0 

one obtains 


= 2p tcRLCq 




lir;2 

Po + ^ t2 


e* 5' dri 
b 0,0 b 0,l 


( 5 . 26 ) 


and 


W B = 2pnRLC 


$ 


h „(o) 


J °JL° 


0 ! 

2 R< 


£t 2 


(5.27) 


The function 


-0,2a^ 


is explained at the end of the appendix. We 


did not numerically determine the coefficients at t 2 in equations ( 5 . 26 ) 

and ( 5 . 27 ). Because of Oq = i ^o'q( 0)(= l/ /fn = 0.5642) the law stated 

above W = W - W is valid also in this limiting case n = 0 of 
R D D„ 


the sudden start of motion, for times immediately after the start of 
motion. Here in particular Wp =0 (d’Alembert's paradox). 
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APPENDIX 


A FEW CALCULATIONS REGARDING THE DEVELOPMENT OF THE BOUNDARY LAYERS 


A few calculations will be given in this appendix which yield, 
among other data, those required for the preceding investigation 
(section 5 ) concerning the basic functions £ n -^(q) of the unsteady 
boundary layers; for the rest, they represent Merely an extension of 
the related calculations by Blasius. We give these calculations apart 
from the previous considerations, first, because they would have 
disrupted the connection there, and second, because the datu and tables 
attained are of interest in their own right. 

As assumed in section 5, let a velocity proportional to t n (n = 0) 
be im parted to a body from a state of rest relative to the surrounding 
indefinitely extended fluid, beginning at t = 0. The potential- 
theoretical circumferential velocity U(x,0,t) = Uo(x,t) then has the 
form (5.1). For calculation of the boundary layer development from 
t = 0, if a generalization is made of the series set up by Blasius 
for n = 0 and n = 1, the expression 


t(x,T!,t) = 2^Tt Y~ t n+ ^( n+1 ) X (x,q) 
with (l) 

T) = y/2 fjvt 


for the stream function of the boundary-layer flow is obtained and 
one obtains for the functions X n ^ by substituting equation (l) 

into the boundary layer differential equation (l.l) a system of 
differential equations solved by recursion; we limit ourselves here to 
the two first equations of this system which read 


83x 


n,0 


S 2 * 


3tt 




n,i 




+ 2q 


+ 271 


n,0 


8 2 x 


^1 

n,! 


br\ 


8 X n 0 

bn — = - 4nq(x) , 

8ti 

8 x /8 x _ 8 c x _ 

■ / x n,l 1 I n,0 n,0 

4(2n + 1) ■ - - = k 1 ' * 


8t] \ 8tj 8x87] 

a \,0 ^Xn.O 


>( 2 ) 


8x 8 t]‘ 


- q 4' • 
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We assume first n (later 2n) to be an integer. Then the . solutions 
may be represented with the aid of Hermite polynomials. With the 
statements 


X 

n,0 


<l(x)£ (ti), 

n,0 


> 


X n,l = 


(3) 


one obtains from equation (2) the ordinary differential equations 


+ aiC* * 

n,0 n,0 


t ' ' ' + 2nt* ' 
b n,l 1 n,l 


- 4n£* = - 4n. 

n,0 


(V 


- M2n + l)t; (1 = lv(C,o 2 - £n, 0 5n,0- l)- (5) 


The boundary conditions to be satisfied by t _ and £ , 

^n^O ^n,l 

formulated in equation (5-*0 


are 
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In the case of a plane wall moving in its own plane. 


= 2 qt n £ n qC 1 !) represents the complete solution .(because 

of q = const . ) , not only in the boundary-layer theory approximation 
set up here, but in strict fulfillment of the complete Navier-Stokes 
" 1 20 

equations. The calculation of (; may take place as follows. The 

n *° 


temporary splitting-off 


= 1 " e ” T1 ■ < P n ( T l) 
n,0 11 


( 6 ) 


transforms equation (4) into 


'Pn " ^n “ 2 ^ 2n + 1 ^ ( P n = 0 


( 7 ) 


20-r 


Because of U Q = U Q (t), u = u(y,t) and v = 0, the Navier- 
Stokes equations are for these motions simplified to 


Su 

S 2 u 

SlL 

= V — 5 

+ — 0 - 

St 

dy 

St 

u 

and the 

corresponding solutions of the 


problem of heat conduction has been known for a long time; it offered 
one of the few possibilities of attaining exact solutions of the 
Navier-Stokes equations. For the rest, one can see for the present 
problem that the first term qt£^ q of the series development 

following from equation (l) for u as a solution of the above equation 
approximately satisfies the boundary-layer equation in the sense that 
only the terms of highest order are taken into consideration for sma ll 
times after start of the motion, whereas the quadratic inertia terms 
are neglected. The iterative improvement of this first approximation 
for small times then yields step by step the ascending terms of the 
series we set up formally at the outset. This consideration led 
Blasius, (elsewhere), at the time to his special series formulations 
for n = 0 and n = 1. 
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For every integral 2n ^ 0 the general solution of this differential 
PI 

equation is 



(Ci, C 2 are integration constants). As is well known, the Hermite 
polynomials H^x) are given by 


H (x) 
m v ' 



°< k < £■ 


<- »*{&)*&■ 


(2x) 


m-2ic 


( 9 ) 


(in the form originally given by Hermite). For further use we also put 
H m (ix) = i“8^(x), thus 


H (x) 
m 




d 

Vdxy 


0 < 


& 


2 k): 
k ! 


(2x) 


m-2K 


( 10 ) 


furthermore 


® (t) 


0 



V x) 



( 11 ) 


^Compare, for instance, E. Kamke "Differentialgleichungen: 

Lo sung smetho den und Losungen, I Gewohnliche Differentialgleichungen" 
(Differential equations: Methods of solutions and solutions, I. 

ordinary differential eauations) Leipzig 19^2, Part C, No. 2,4l, and 
put there x = i ^2 q . 
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(so that in particular 


1+0 = 0 
0 


represents the error integral, 0 ^ 


the error function). By using these expressions, a simple recalculation 
from equation (8) gives the general solution of equation (4) as 


V Oi) 

n,0 


1 - H (n)|c + C*o (tj) + 

2n j_! 2 0 _ 


C*®^) 



(- D K f 


2n-K 


(r 1 )H K _i( T i ) 


( 12 ) 


with C 2 * = Cg. 

= 0 because of £ ' (<») = 1 and C 2 * = -n.'/(2n).' because 

of £' (0) = 0 and H2n(0) = (2n)i/n.* . ["The polynomial sum in 

n l* 

equation (12) disappears for T) = 0 and integral n term by term, 

since either 2n - k or k - 1 is an odd number! . Thus the desired 

PP ~ 

particular integral of (4) reads 


5 ' (n) = 1 + 7-=?- 

n,0 (2n) 


2n 


H 2n (Tl)* 0 (n) - SqOl) XT ( 2n )(- 1 ) K H2n- K ( T l) H K -l( 1 l) 

K =1 ' K ' 

(13) 


Because of $> (tj) = (- l) K_ ’ 1 ~0- 1 (q)H , (t)) the solution may also be 

written in the following form which is more elegant than equation ( 13) 
but less serviceable for practical calculation 


h n<1> - 1 


n: 


"n,0 


(2n)i it=0 




(13a) 


22 Th 


e direct and elementary derivation of the solution for 


integral n 3 0, and therewith for integral 2n 3 0, compare below, 
fails to work if n does not have this property. But in that case, 
too, the solutions are easily found if one makes use of the analogy to 
the corresponding solutions of linear heat conduction (compare 
footnote 20 to this report) and represents the solution according to 
the singularity method. 
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According to equation (13) one has in particular 


• = 1 + <c 0 (,l) = 

S' Oi) = 1 + (2^1 2 + i)® 0 (ii) + n^Oi), 

^ 0 (n) = 1 + 3 C (4t)4 + 12r[2 + 3 ) ®o (n) + (2T]2 + 

^3 0 ^ = 1 + ^ + + 90t \ 2 + 15)® 0 (tj) + > 

(W*+ 28t ] 2 + 33)TJ« 1 (r|’)j , 

otn) = 1 + ^-^|i6t|8 + 224t^ ^ + 840ti^+ 840tj2 + 105) ^(b) + 
(8t] 6 + 108t 1 4 + 370n 2 + 279)11^(11)]. 


(13h) 


Thence one obtains by elementary quadratures t ^Likewise 

expressed by 4> Q and ^ with polynomial coefficients. The numerical 

evaluation is reproduced in table 1 2 3 , Because of the special impor- 
tance of the solutions £’ (ri) as boundary- layer profiles u/qt n 

n ,0 

in the case of the plane wall (compare above) we represented them in 
figure 3 . [Owing to 6 * - 2 f \[yt ^ljm_ (t| - 0 (q)) = 2a n tbe 3 ' tream - 

line displacement thickness 5 * can easily lae taken from the numerical 
cal table 1 ; compare -also table 2.1 


Since it follows directly from equation ( 7 ), by single differ- 
entiation, that the general integral cp -> is the first derivative of 

n+ 2 

the general integral <P n with respect to rj, it is, with the 
boundary conditions satisfied, easy to find as expressions for the basic 
functions £ (t) ) with the integral n 

n+i -,0 


'4° ‘ 1 ' C,o< 0) P’° + ^ 


(l4) 
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All numerical calculations were performed by Miss Ursula Ludewig. 
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Therein 


5 " ( 0 ) = 2 2 * 1 
n,0 (2n)j a/ji 


(15) 


In order to calculate the fundamental f met ions of the first order 
^n the total course of ^(t]) mus t he known (according to 

equation (5))- Only the additional knowledge of (0) is required 

n jl 

for the problem which is of foremost interest, the question regarding 
location and time of a possible separation of the boundary layer in 
first approximation. 


For because of 





(0) + qq't n+1 ^ 1 (0) + 




one obtains in first approximation the connection 



(16) 


for location x and time t of the separation. On the other hand, 

we are interested in £ ( 00 ) (compare section 5), with a view to the 

n,l 

calculation of the displacement thickness. These two data can be 
determined without solution of the differential equation (5) by a 
well-known method as follows: $(t|) is assumed to be a function of t} 

in 0 = Tj = 00 provided with the continuity properties required for the 
following calculation. By partial integration one obtains the following 
relations. If L n and M n are the differential operators 


•,2 

L n = -^p + 2t ) — - M2n + 1), 

dT] 
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one has 



We choose, therefore, -d = -3 (t)) so that 

n,0 


m Pa = 0 

B L n '°J 

with (17a) 


■a (0) = - 1, a («) = 0 

n,0 n,0 


If then equation ( 5 ) and the boundary conditions valid for 
are taken into consideration the result is 



(n) 



On the other hand we choose i3 = (ri) so that 

car 1 # 


■ 1 

with 

$ (0) = 0, $ («,) finite 

n^«> n^oo 


A 


(17b) 


(18a) 
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Then one obtains, in analogy 


£ (») = [ t 3 (1)4/^* 2 f f" 

n,l JO n,oo 4 ^n,0 “ %0h,0 " i/ dT l. 


(18b) 


It is easily confirmed that 


n , 00 


(n) = 


-1 


2(4n + 3) 




(19) 


furthermore by comparison with equation ( 7 ) 


v ( ti) = 


p2n+l,0^ " X ] J 


( 20 ) 


Therewith the desired functions $ n 

n,0 

the known basic functions of zero order 


and 

i 


2tj+1 ,0 . 


are traced back to 


Numerical evaluation yielded for n = 0, 1, 2, 3, and 4 all 
together the data here of interest given in table 2. It also shows the 
numerical values of a n . For the n-values 1, 2, 3, and 4 one finds 
the law (5.24) confirmed. A general proof of this law may be produced 
with a few calculations on the basis of the known expressions 
for £ n o(h)* ^ mucd simpler proof of the relation (5.24) will be 

presented below. As mentioned above, the expression 



is the strict solution for the boundary- layer profile on a plane wall 
with U = qt n (q. = const.) outside of the boundary layer. The wall 
shearing stress according to the momentum theorem of the boundary 
layers is generally 


U 


8x ^0 


pudy - — 
ox 


pu 2 dy ■- 


Su , „ dp 

p _ dy - 5 

3t ox 


( 21 ) 
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thus for the above flow, due to the velocity distribution u(t) ,t) being 

independent of x as well as due to = - p and with 

ox dt 


- ^oO 1 - 4,o (,l) ] 


T n = p r (u - u)dy = p ^-(U5*) = pq(2n + l)a Jv 2. (22) 

u dt JO St n r 


On the other hand 


T _ 


= M- 


Su(0,t) 
Sy 


Vv + n "| 


p ^,0 (0) -¥ * 


(23) 


and thus in combination with equation ( 15 ) , as asserted, 

22n+l/ n t\2 

2(2n + 1 ) 0 ,, = ^ 0 (°) = (2n).* ^ ’ 


( 5 . 24 ) 


In the case of the flow about a circular cylinder moved rectilinear ly 
out of a state of rest, investigated in section 5 [compare equation (5.8)^] 
one has q(x) = 2Cq sin x/R and therefore -q'(x) max = 2Cq/R 

for x/R = jt. Thus the separation starts according to eqmtion (l6) 
at the rearward stagnation point at the time 




1 

n+1 


(compare table 3). The distance covered by the cylinder during that 
time is 


S A = C O t A n+1 /( n + a ) = R £" (0)/2(n + lK" (0) 

A yJ A / n,0 n,l 
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Finally we want to give a few indications where to find further 
data regarding the basic functions of the plane nonsteady boundary-layer 
flow. The series development (5.2) written down up to the third term 
reads 


* - 2 4 * + tn+1< l ,E n,l + t2(n+1, ( q ' 2£ n,2a + 91 " E n,2b) + "D* 


Blasius gives, in addition to the functions £q (tj) and q( 1 i) 
calculated above, the rigorous solutions £ (tj) and ^ 

Beyond that, he calculates the numerical values £" g (0) and £" g, (0) 

which are of interest for the determination of the separation. 

S. Goldstein and L. Rosenhead*^ give the exact expressions for ^o,2a( T l^ 
arid ^ These integrals were, by the way, numerically determined 

before by Boltze^ on the occasion of treatment of the corresponding 
problem n = 0 of rotationally symmetrical flows which seems to have 
escaped the attention of the authors. 


Translated by Mary L. Mahler 
Rational Advisory Committee 
for Aeronautics 


^Compare footnote 6 of this report. 


''Compare footnote 5 of this report. 



TABLE 1 


■ 

£0,0 

£0,0 

£5,o 

£1,0 

? i,o 

t" 

%0 

£2,0 

^2,0 

t" 

^2,0 

h,0 

? 3,0 

3,0 

5 4,0 

S 4,0 

t" 

M 

0 

0.1 

0.2 

0.3 

o.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

1.1 

1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

1.8 

1.9 
2.0 
2.1 
2.2 

2.3 

2.4 

2.5 

2.6 

2.7 

2.8 

2.9 

3.0 

3.1 

3.2 


0 

0.1125 

0.2227 

0.3287 

0.4284 

0.5205 

0.6039 

0.6778 

0.7421 

0.7969 

0.8427 

0.8802 

0.9103 

0.9340 

0.9523 

0.9661 

0.9763 

0.9838 

0.9891 

0.9928 

0.9953 

0.9970 

0.9981 

0.9989 

0.9993 

0.9996 

0.9998 

0.9999 

0.9999 

1.0000 

1.0000 

1.0000 

1.0000 

1.1284 

1.1172 

1.0841 

1.0313 

0.9615 

0.8788 

0.7872 

0.6913 

0.5950 

0.5020 

0.4151 

0.3365 

0.2673 

0.2082 

0.1589 

0.1189 

0.0872 

0.0627 

0.0442 

0.0305 

0.0207 

0.0137 

0.0089 

0.0057 

0.0036 

0.0022 

0.0013 

0.0008 

0.0004 

0.0003 

0.0001 

0.0001 

0.0000 

Hi 

HSR » J 

0 

0.20 64 
0.3773 
0.5172 
0.6301 
0.7201 
0.7910 
0.8459 
0.8880 
0.9197 
0.9432 
0.9604 
0.9728 

0.9816 

0.9878 

0.9920 

0.9948 

0.9967 

0.9979 

0.9987 

0.9992 

0.9995 

0.9997 

0.9998 

0.9999 

1.0000 

1.0000 

1.0000 

1.0000 

2.2568 
1.8793 

1.5464 

1.2569 
1.0085 
0.7986 
0.6237 

0.4804 

0.3647 

0.2728 

0.2010 

0.1459 

0.1042 

0.0733 

0.0507 

0.0345 

0.0231 

0.0152 

0.0098 

0.0063 

0.0039 

0.0024 

0.0015 
0.0008 
O’. 000 5 
0.0003 
0.0002 
0.0001 
0.0000 

0 

0.0138 

0.0506 

0.1049 

0.1721 

0.2489 

0.3327 

0.4215 

0.5138 

0.6087 

0.7052 

0.8030 

0.9015 

1.0006 

1.1000 

1.1996 

1.2994 

1.3993 

1.4992 

1.5992 

1.6991 

1.7991 

1.8991 

1.9991 

2.0991 

2.1991 

2.2991 

0 

0.2638 
0.4638 
0.6139 
0.7251 
0.8066 
0 . 8656 
0.9077 
0.9374 
0.9581 
0.9724 
0.9820 
0.9884 
0.9927 
0.9954 
0.9972 
0.9983 
0.9990 
0.9994 
0.9996 
0.9998 
0.9999 
0.9999 
1.0000 
1.0000 
1.0000 
1.0000 

3.0090 

2.2941 

1.7298 

1.2896 

0.9501 

0.6916 

0.4972 

0.3529 

0.2473 

0.1709 

0.1166 

0.0784 

0.0520 

0.0340 

0.0219 

0.0139 

0.0087 

0.0053 

0.0032 

0.0019 

0.0011 

0.0007 

0.0004 

0.0002 

0.0001 

0.0001 

0.0000 

0 

0.0162 

0.0584 

0.1189 

0.1923 

0.2745 

0.3628 

0.4552 

0.5502 

0.6471 

0.7451 

0.8439 

0.9432 

1.0427 

1.1425 

1.2423 

1.3422 

1.4422 

1.5421 

1.6421 

1.7421 

1.8421 

1.9421 

2.0421 

2.1421 

0 

0.3067 

0.5245 

0.6718 

0.7835 

0.8564 

0.9059 

0.9390 

0.9610 

0.9753 

0.9846 

0.9905 

0.9942 

0.9966 

0.9980 

0.9988 

0.9993 

0.9996 

0.9998 

0.9999 

0.9999 

1.0000 

1.0000 

1.0000 

1.0000 

3.6108 

2.5762 

1.8184 

1.2695 

0.8763 

0.5978 

0.4031 

0.2685 

0.1766 

0.1147 

0.0736 

0.0466 

0.0291 

0.0179 

0.0109 

0.0065 

0.0038 

0.0023 

0.0013 

0.0007 

o.ooo4 

0.0002 

0.0001 

0.0001 

0.0000 

0 

0.0182 

■0.0646 

0.1298 

0.2075 

0.2934 

0.3845 

0.4790 

0.5757 

0.6736 

0.7724 

0.8717 

0.9713 

1.0710 

1.1709 

1.2708 

1.3708 

1.4708 

1.5708 

1.6707 

1.7707 

1.8707 

1.9707 

0 

0.3415 

0.5710 

0.7235 

0.8237 

0.8889 

0.9307 

0.9573 

0.9740 

0.9844 

0.9907 

0.9946 

0.9969 

0.9982 

0.9990 

0.9994 

0.9997 

0.9998 

0.9999 

0.9999 

1.0000 

1.0000 

1.0000 

4.1266 

2.7858 

1.8608 

1.2296 

0.8035 

0.5192 

0.3316 

0.2093 

0.1305 

0.08o4 

0.0489 

0.0294 

0.0175 

0.0102 

0.0059 

0.0034 

0.0019 

0.0010 

0.0005 

0.0003 

0.0002 

0.0001 

0.0000 
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TABLE 2 


n 

“n 

Sn,0«» 

Sn,lH 

- £S,!(0) 

0 

0.564-2 

1.1284 

0.418 

1.607 

1 

0.3761 

2.2568 

0.138 

0.963 

2 

0.3009 

3.0090 

0.072 

0.756 

3 

0.2579 

3.6108 

0.046 

0.632 

4 

0.2293 

4.1266 

0.033 

0.552 


TABLE 3 

START OF SEPARATION t A AND DISTANCE TRAVELLED S A IN THE 
CASE OF A CIRCULAR CYLINDER MOVED RECTILINEARLY 
OUT OF A STATE OF REST 


n 

1 

Cp 1 n+1 

tiM 
A l R J 

S _A 

R 

0 

0.351 

0.351 

1 

1.082 

0.586 

2 

1.258 

O.663 

3 

1.300 

. 0.714 

4 

1.302 

0.748 
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Figure 2.- Pressure drag W £>, friction drag Wr, and total drag 

W = W-q + W R of the circular cylinder for very small times after start 
of the motion for uniform acceleration out of a state of rest. 






